Entanglement in spin chains and lattices with long-range Ising-type interactions 
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We consider N initially disentangled spins, embedded in a ring or d-dimensional lattice of arbitrary geometry, 
which interact via some long-range Ising-type interaction. We investigate relations between entanglement 
properties of the resulting states and the distance dependence of the interaction in the limit N — * oo. We provide 
a sufficient condition when bipartite entanglement between blocks of L neighboring spins and the remaining 
system saturates, and determine Sl analytically for special configurations. We find an unbounded increase of 
Sl as well as diverging correlation and entanglement length under certain circumstances. For arbitrarily large 
N, we can efficiently calculate all quantities associated with reduced density operators of up to ten particles. 
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The investigation of entanglement properties of strongly in- 
teracting many body systems has proven to be a fruitful ap- 
proach. Entanglement was shown to indicate quantum phase 
transitions lUlSIIlBl an d the presence of long-range correla- 
tions even in systems with gapped Hamiltonians |5|]. In den- 
sity matrix renormalization group (DMRG), a powerful nu- 
merical method which allows one to treat spin chains of up to a 
few hundred particles with high accuracy, the investigation of 
the role of entanglement has allowed one to understand 001 
and overcome |01 limitations of the method. Standard DMRG 
can treat spin systems with a bounded amount of bipartite en- 
tanglement, measured by the entropy of entanglement Sl be- 
tween blocks of L neighboring spins and the remaining sys- 
tems (e.g. spin chains with short range interactions). The 
generalized method of Ref. j^t can also handle spin systems 
arranged on two (and higher) dimensional lattices, where Sl 
scales with the surface of the block, Sl oc L 1 / 2 . 

Despite these significant developments, many spin sys- 
tems in two and three dimensional setups remain untractable, 
among them disordered systems or systems with long range 
interactions where Sl oc L. Such spin lattice systems with 
long-range interactions occur naturally in several quantum 
optical setups. There, Ising-type interactions are induced by 
other interactions with a characteristic distance dependence. 
Examples are the internal states of neutral atoms in an optical 
lattice that interact via an induced dipole interaction |9|] (see 
also Jioll ). or ions stored in microtraps where interactions are 
induced by pushing the ions dependent on their internal state 
such that they feel a different Coulomb potential fTIHT2lfl3ll . 
In this letter, we analyze spin chains and spin lattices in ar- 
bitrary dimensions with long-range Ising-type interactions. 
Despite the failure of known methods, the restriction to Ising- 
type interactions allows us to analyze both static and dynami- 
cal entanglement properties of the system in great detail. We 
find relations between entanglement properties of states |^ t ) 
— resulting from evolution of the system initially prepared in 
some product state l^o) under the Hamiltonian H for time 
t — and the distance dependence of the interaction. We obtain 
information about the dynamics of entanglement and, at the 
same time, a large family of multipartite entangled states with 
rich entanglement properties. 



Our results are based on a description of the states |\I> t ) in 
terms of generalized Valence Bond Solids (VBS) [14]. In this 
picture, we can efficiently calculate the reduced density op- 
erators of a small number L < 10 of arbitrary spins, even if 
the total number of spins N is large (a standard PC can eas- 
ily handle N = 10 5 particles). We can hence determine all 
quantities associated with reduced density operators of small 
subsystems, including e.g. higher order correlation functions 
or (bounds on) bipartite entanglement Sl- We emphasize that 
for general pure states, the calculation of reduced density op- 
erators is a highly non-trivial task due to the exponential scal- 
ing with the system size 7Y. For certain distance laws, we can 
describe the scaling of block-wise entanglement with the size 
of the block L in the limit N — ► oo. In one dimension, we 
obtain a sufficient condition when block-wise entanglement 
saturates, which is the case whenever the distance dependence 
of the interaction strength scales as 1/ r a with a > 1. For spe- 
cial configurations, we calculate Sl analytically for all L and 
show that Sl can indeed grow unboundedly and proportional 
with L. This is in contrast to en tang lement properties of 1D- 
VBS states recently analyzed in lfl5ll . where Sl is bounded by 
two. Finally we find that the correlation length diverges under 
certain circumstances, even if Sl saturates. 

In our model, we consider N spin 1/2 systems (qubits) with 
pairwise interactions, described by an Ising-type Hamiltonian 
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We assume that the spins are arranged on a d-dimensional 
lattice with fixed geometry and are initially completely po- 
larized in x-direction, i.e. |^o) = |+)® > where |+) = 
l/\/2(|0) + |1)). The methods we develop can also describe 
disordered systems with random coefficients f(k, I) and can 
take arbitrary (product) input states into account. We are in- 
terested in (entanglement) properties of the state 
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We consider the situation where the coupling between spins 
obeys a certain distance law, in the sense that the coefficients 
f(k,l), describing the strength of the coupling, only depend 
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on the distance rki = ||A; — Z|| between particles k and I, 
.f(k, l) = f{Tki)- In the example of ions stored in microtraps 
lflll. fl3l1 one finds for instance f{ru) = r kl 3 113ll - 

Description in terms of Valence Bond Solids: The uni- 
tary evolution operator U (t) = e~ ltH in (0 can equivalently 
be described by a product U = Y\ k l Uki of commuting con- 
trolled phase gates E/jy = diag(l, 1,1, e l< ^ kl )ki acting on pairs 
of qubits, where <pki = f(rki)t- If U acts on a completely po- 
larized state |\P ) = 1+)®^ and 4> k i £ {0,7r}, the resulting 
states are graph states flfil fl7ll . Hence, we will refer to states 
with arbitrary ipki (produced from |+)® ) as weighted graph 
states. Verstraete and Cirac recently proposed a description 
of graph states in terms VBS Q. Instead of maximally en- 
tangled pairs of qubits, we use pairs where the degree of en- 
tanglement depends on (p^. In this generalized VBS-picture, 
we develop a description of all weighted graph states and ex- 
tend it to all states produced by the interaction Hamiltonian H 
acting on any initially unentangled state. 

Each qubit k of a physical state \ty t ) is replaced by N — 1 
(virtual) qubits k\ . . . fcjv-i- The VBS state \^> t ) with corre- 
sponding HUbert space Ti = [(fT 2 ) N ~ 1 ] N is given by a ten- 
sor product of N(N — l)/2 independent, non-maximally en- 
tangled pairs of qubits \xkilj) = Ukd\+)k l \+)i 3 shared be- 
tween virtual qubits of parties k,l, where each ki,lj 
appears only once. Up to a normalization factor, we ob- 
tain the corresponding weighted graph state \ty t ) = U\+)® N 
from |^ t ) = ®IXfei?i) by performing local projections 
Pk = |0fe)(0?| + |lfc)(ld onto two dimensional subspaces 
at all locations k, where k = k% . . . fcjv-i and M) = 
\mm . . . m). The VBS-like state \^>t), together with the pro- 
jection <^) k Pk, thus provides an equivalent description of the 
state \^t)- We can generalize this description to arbitrary 
product input states \ip\ , . .tpjsf). In this case, the (unnormal- 
ized) VBS-like state is of the form \^ t ) = ® fc ; \xk z ij) with 
|XWi) = Uki\ "-jftpktkil N -^I)i 3 , where Iv?) = a|0)+j9|l) 
and| irfip) = \fa\0) + v^3|l). In the following, we use this 
description to determine reduced density operators pa- For 
the sake of simplicity, we consider weighted graph states, i.e. 
states arising from input states 1+)®^. 

We denote by A an arbitrary subset of the N qubits, and 
we call the set of remaining qubits B. Because all Uki com- 
mute and because unitaries in B do not influence pa = 
trsd^t) (^t\), we can write 

Pa = II U u tvB\*' t )(V t \Ul v (3) 

k,lGA 

with |^) = Y\ keA ,ieB u ki\^o)- We now determine p' A = 
tr B |*i)(*i| in the VBS-picture, i.e. we start with the state 

l*J> = (g) Uki\+) kl \+) h . (4) 
keA,ieB 

For the following argumentation it is crucial that Uki = 1 for 
k, I € B as is the case for \^f' t ) but not for |$t). Performing 
the projections Pi on all particles I e B (but not in A) leaves 
us with a state of the form (S^esIIIfceA Uki\+) k \+)i\- We 



have a tensor product of \B\ states where in the I th state par- 
ticle I € B is entangled with a virtual system A. For each of 
these states we can independently calculate the reduced den- 
sity operator with respect to A. We simply trace out particle I 
and obtain 

p'aQ) = U\+)a(+\ + \xi}a(xi\), (5) 

with \ X i)a = <8> feeA (l°> + e-^'\l))/V2. Now, we per- 
form the projections Pk for all k <G A. The resulting density 
operator p' A is (up to normalization) given by the Hadamard 
product of all density operators p'a(1), where the Hadamard 
product of two matrices corresponds to component-wise mul- 
tiplication in the computational basis. The reduced density 
operator pa is then obtained from p' A by taking into account 

interactions within A leading to pa = life ieA ^uPa^m- 
nally, we must normalize the resulting state. 

Computable quantities: The method outlined above pro- 
vides an efficient way to calculate reduced density operators 
for all weighted graph states. The computation time is linear 
in the number \B\ of particles in the remaining system (but 
exponential in \A\), as opposed to an exponential scaling in 
N = \A\ + \B\ of computation time and memory cost for 
general pure states. Hence, for arbitrary large systems, all 
quantities that depend on the reduced density operator of a 
small number of qubits can be calculated efficiently. For in- 
stance, from pa of one and two qubits, we can determine all 
two-point (and also higher order) correlation functions Q a ' „, 

Q& = (°P*P) (6) 

lower and upper bounds on the localizable entanglement El 
0, the entanglement of formation between pairs of parti- 
cles, as well as the multipartite entanglement measure -Emw 
lTl8ll . The maximal classical correlation between two 

particles is given by the largest singular value of the ma- 
trix 0. The localizable entanglement E k L ,x is the 
maximum amount of entanglement that can be established 
between a pair of particles k, I, on average, by perform- 
ing local measurements on all other particles. The relation 
QmL < e l 1 < ^^holds 0, where E k / is the con- 
currence of assistance Il9ll . The measure Emw is given by 
Emw = 2[1 - l/NEfc tr 0°I)] H- In addition, the (bipar- 
tite) entanglement between blocks of a small number L < 10 
of neighboring spins is measured by the entropy of entangle- 
ment Sl, that is the von Neumann entropy of the reduced den- 
sity operator pL, Sl = tr(/j£ log 2 pl)- Clearly, < Sl < L, 
where Sl = L indicates maximal entanglement between the 
blocks. For blocks larger than 10 qubits, we make use of the 
strong subadditivity of the entropy to derive upper bounds on 
Sl- By breaking a block of size L into n sub-blocks Li of size 
| Li | = L/n, we obtain 

n— 1 n— 1 

S(pl)<5>( PLi ,l i+1 )-E 5 (^)' W 

i=l i=2 

where pL i: L i+1 is a sub-block of length 2L/n. The unitary 
operations Uki ,k,l€A,do not change the entropy S(pa) and 
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hence the reduced density operator p' A can be used directly 
since S(pa) = S(p' A ). Nevertheless, upper bounds on the 
entropy are different for pa and p' A , where the latter turn out 
to be more stringent and will hence be used in the following. 

The fact that the total reduced density operator is given by 
the Hadamard product of reduced density operators with re- 
spect to all particles in the system B can be exploited to prove 
monotonicity properties of S(p' A ). For a fixed size \A\, we 
add one particle j to B. The reduced density operator is up- 
dated by Hadamard multiplication with p' A (j). From theo- 
rem 5.5.12 in 12011 follows that the eigenvalues of the result- 
ing density operator are majorized by the eigenvalues of the 
initial one, which implies that the entropy increases l2lll . As 
a consequence we obtain lower bounds on the entropy of en- 
tanglement S{p' A ) when we take into account only a subset 
BcBof all particles (and ignore the other particles in B). 

Static properties of resulting states: We apply these re- 
sults to determine (static) entanglement properties of the state 
\^t) for some fixed time t < tt and for different distance 
laws f(rki) = r ki a j a > 0- Figure [Ha) shows the maximal 
two-point correlation Q\^ ayi in a chain of N = 10 5 particles 
as a function of of the distance between particles ||i — 
We observe that correlations decay slower than exponential. 
Therefore, the correlation length £ and also the entanglement 
length £s diverge 01 ■ This indicates long-range quantum cor- 
relations for all power laws, as we find that only exponential 
fall-off functions f(k, I) = e~ Krkl lead to a finite correlation 
length. Figure Qb) shows the scaling of the entropy of en- 
tanglement Sl with the block size L for different power laws. 
Exact values are plotted for L < 10, while upper bounds (cor- 
responding to \Li\ = 4) are plotted for L > 10. The up- 
per bound on Sl seems to grow unboundedly for a < 1/2, 
whereas Sl saturates for a > 1. For a > 1 the system thus 
contains a bounded amount of entanglement Sl, but has a di- 
verging correlation length £. 
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FIG. 1: Spin chain with N = 10 s , t = 0.3tt and f(r ki ) = r~ a for 
different a. (a) Maximal two point correlation Q^ax as a function 
of the distance \\i — j\\ for a = 1/3 (dashed) and a — 3 (solid), 
(b) Exact values [upper bounds] of entropy of entanglement Sl as 
function of block size L for a — 1/3 (A [o]) and a = 3 (op]). 

The saturation of Sl for a > 1 can be proven analyti- 
cally when we take the limit of an infinite chain (N — > oo) 
and afterwards let the block size also go to infinity. Both 
steps involve infinite products and sums, and we must con- 
cern ourselves with convergence/divergence of these products 
and sums. In an infinite chain, the reduced density opera- 
tor for one qubit is p A = ±[1 + c|0)(l| + c*|0)(l|] with 



c = n£Li cos ^6-^=1 <?W 2 an d (j) lk = r -«t IH. The 
eigenvalues of p A are given by (1 ± |c|)/2, so we can omit the 
phase of c. We write YlkLi cos^ik = e^-^ 1 lncos <^ifc. Taylor 
series expansions and Cauchy's integral criterion tell us that 
the sum in the exponent converges for a > 1/2 and diverges 
to — oo for < a < 1/2. Hence, the infinite product is non- 
zero in the first case and zero in the second case. The entropy 
of entanglement Sa is thus smaller than one for a > 1/2 
and equals one for smaller a. In other words, for slow fall- 
off functions (strong long-range interactions), the entropy of 
a single particle is maximal, independent of the time t > 0, as 
infinitely remote regions still influence the qubit we consider. 
To take the limit L — > oo in the case a > 1/2 we use a 

bound given by S(p L ) = S(p' L ) < Y^Li ^(p'lJ for sub " 
block sizes \Li\ = 1. In this case, the convergence proper- 
ties of the infinite sum of single particle entropies (each itself 
given by an infinite product) can be determined by using again 
Taylor series expansion and Cauchy's integral criterion. In the 
limit N —> oo and L — > oo, the upper bound converges to a 
constant value for a > 1 and hence the exact value of the 
entropy Sl also saturates as a function of L for power laws 
f( r ki) = rZ a with a > 1. We can generalize this result to d— 
dimensional lattices. When considering blocks of L particles 
contained in a d-dimensional ball, Sl can at most grow like 
the volume of that ball, whereas we find that for a > d the 
upper bound on Sl grows at most like the surface of the ball. 
The proof is similar to the one-dimensional case. 

For special cases, we get a complete analytic description of 
the entanglement properties. In the following, we consider an 
interaction Hamiltonian with a fixed interaction length A and 
constant interaction strength, i.e. f{rk{) = 1 if rjy < A and 
zero otherwise. For t = tt, the resulting states are spe- 
cial instances of graph states flrjll . We denote a d-dimensional 
quadratic block of size L = a d neighboring spins by A and 
the remaining system by B. We measure the bipartite entan- 
glement between A and B by the entropy of entanglement Sl ■ 
For graph states, Sl is given by the binary rank of the adja- 
cency matrix Tab between the quadratic block A and the rest 
B 11711 . If the lattice is large enough to contain not only the 
block A but also the larger block of size (a+2A) d with A in its 
center, then no boundary effects have to be taken into account. 
We can inductively show that the matrix Tab has maximal 
rank by considering different layers Ak with geometric dis- 
tance k to B, starting with k = 1. Hence, Sl is simply given 
by the number of vertices within A that are connected to the 
rest B, so Sl = a d — [a — min(2A, a)] d with a = \J~L. In 
a general situation, no such simple rule to calculate Sl holds 
(in contrast to what is suggested in Q). A counterexample 
is given by a state with TV = 4, ^13 = (^ 14 = <^ 2 3 = <^24 = tt, 
which has S{fi\^) = 1. For an (infinite) chain of particles we 
obtain Sl = min(2A, L). As L increases, Sl saturates at the 
value 2A for any fixed interaction length A. Only if A itself 
goes to infinity as N — > 00, Sl can grow (unboundedly) with 
L when the ratio of the interaction length to the total number 
of particles is kept constant. For a given interaction length 
A, and for any ball containing L particles, only those parti- 
cles that have connections to the remaining system contribute 
to Sl- The entropy Sl equals L if the radius of the hyper- 
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sphere is smaller than A. Otherwise, Sl scales essentially 
like the volume of a surface shell with thickness A , that is 
Sl oc AL' <i ~ 1 - ) / d . Two-point correlation functions Q^fg turn 
out to be zero for all pairs of particles (except for the case of 
two neighboring particles at the end of a chain), because their 
reduced density matrix pki is the identity fPHl . Nevertheless, 
\^!ir) is maximally connected flol . which means that a Bell 
state between any pair of particles can be obtained by local 
measurements on the rest of the particles. Hence has 
maximal localizable entanglement between any two particles, 
Ei = 1, and, thus, an infinite entanglement length. Finally, 
Emw = 1 for a H such states since the reduced density opera- 
tor of each single particle is also the identity. 

Dynamics of entanglement: We return to the case of a 
chain of particles and to Hamiltonians with arbitrary distance 
dependence f(k, I) = rZ, a and consider the dynamics of en- 
tanglement, that means the change of entanglement and corre- 
lations of the state \fy t ) with time. The scaling of the entropy 
with block size L is essentially still governed by the the spe- 
cific form of the distance dependence for any finite t, because 
infinitely remote regions still influence a subsystem in a simi- 
lar way as discussed before. 
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FIG. 2: Dynamics of entanglement for chain with N = 10 and 
f( r ki) = r^ 3 . ( a ) Entropy of entanglement Sl for blocks L = 1 
(bottom) up to L — 7 (top), (b) Two-point correlations Q^ax f° r 
||i — _7 1| = 5 (solid) and Bmw (dashed). 

For large times t, more and more of the interaction phases 



4>ki = f(k, l)t start to oscillate (as they are effectively taken 
modulo 7r) and approach in the limit of large t a (quasi)- 
random distribution. In the limit of an infinite chain and 
t — > oo, the entropy of the reduced density operator of any 
finite group A of particles is maximal, S(pa) = \A\. This 
can be seen by considering the off diagonal elements of re- 
duced density operators, which all contain infinite products of 
cosines of (sums of) random angles. All these products tend 
to zero for N — * oo, leading to a maximally mixed state. For 
a chain of N = 10 5 particles with /(V/y) = r^ 3 the time de- 
pendence of the entropy of entanglement for blocks up to size 
7 is plotted in Fig. |2ja), while Fig.|2jb) shows two point cor- 
relation functions Q max and the multipartite measure Emw ■ 

In this paper, we have investigated entanglement properties 
of states generated from product states by long-range Ising- 
type interactions. For an (arbitrary) total number of particles 
N, and using a description in terms of generalized VBS, we 
could efficiently determine the reduced density operators of a 
few (< 10) particles and hence all quantities which are asso- 
ciated with reduced density operators. For different distance 
laws, we investigated the scaling of block-wise entanglement 
and showed that in ID Sl saturates for /(r^) oc rZ, a for all 
a > 1. We also found diverging correlation- and entangle- 
ment lengths for all power laws. Our methods can also be 
applied to disordered systems, such as quantum lattice gases, 
spinglasses, or the semi quantal Boltzmann gas (introduced 
in 11711 ): they can be extended to describe the dynamics of ar- 
bitrary slightly entangled input states (in the sense of small 
Schmidt measure I23l0 under this Ising-type interaction. We 
remark that related studies have recently been performed for 
harmonic lattice systems in Ref. l24ll . 
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